Macroion correlation effects in electrostatic screening and thermodynamics of highly 

charged colloids 



R. Castaneda-Priego 1 ^ L. F. Rojas-Ochoa 1 ' 2 , V. Lobaskin 3 , and J. C. Mixteco-Sanchez 1 
1 Instituto de Fisica, Universidad de Guanajuato, Loma del Bosque 103, 37150 Leon, Mexico 
2 Departamento de Fisica, Cinvestav, Av. Instituto Politecnico Nacional 2508, 07360 Mexico D. F., Mexico and 
3 Physik-Department, Technische Universitat Miinchen, James-Franck-Str. , D-85747 Garching, Germany 

We study macroion correlation effects on the thermodynamics of highly charged colloidal suspen- 
sions using a mean-field theory and primitive model computer simulations. We suggest a simple 
way to include the macroion correlations into the mean-field theory as an extension of the renormal- 
ized jellium model of Trizac and Levin [Phys. Rev. E 69, 031403 (2004)]. The effective screening 
parameters extracted from our mean-field approach are then used in a one-component model with 
macroions interacting via Yukawa-like potential to predict macroion distributions. We find that 
inclusion of macroion correlations leads to a weaker screening and hence smaller effective macroion 
charge and lower osmotic pressure of the colloidal dispersion as compared to other mean-field models. 
This result is supported by a comparison to primitive model simulations and experiments for charged 
macroions in the low-salt regime, where the macroion correlations are expected to be significant. 

PACS numbers: 82.70.-y, 61.20.-p 



I. INTRODUCTION 

Structure and thermodynamics of charged colloidal 
dispersions became a subject of a renewed interest over 
the last decades due to development of experimental and 
theoretical techniques and accumulation of new data in- 
compatible with classical theories 0, 0, ■ Considerable 
theoretical efforts have been invested into an upgrade of 
existing mean-field approaches with an inclusion of addi- 
tional correlation effects such as counterion or macroions 
correlations, which are missing in the classical Poisson- 
Boltzmann theory. 

Whereas in the aqueous dispersions of micrometer- 
sized particles the correlations of monovalent counterions 
are usually negligible, it is known that they might become 
important for small strongly charged macroions 0, 
The role of counterion correlations has been extensively 
studied by various means, including integral equation 
theories and molecular simulations starting from eight- 
ies 0, IE 0, H E3 El E3 and is currently well understood 
UEl 001 In contrast, macroion correlation effects 
are usually not included in the mean-field approaches and 
therefore not quantified. One can expect these effects to 
be significant in systems with thick double layers, such 
as deionized colloidal dispersions. 

In order to specify our interest in macroion correla- 
tion effects we would like to start from a simple en- 
ergy argument. Charged colloidal suspensions are com- 
posed of a large number of particles of different types. 
If the molecular details of the solvent and dielectric 
discontinuities are neglected, one arrives to the prim- 
itive electrolyte model. On this level, a charged col- 
loidal dispersion is an asymmetric electrolyte consist- 
ing of strongly charged macroions and small countcri- 



*Electronic address: ramoncp@fisica.ugto.mx 



ons. In addition, at least two different species of salt 
ions are usually present. A straightforward application 
of the Debye-Huckel-like mean-field description is usu- 
ally not successful due to strong spatial correlations of 
different ionic species. To deal with the correlations, 
one can attempt to construct an hierarchy of interac- 
tions from a quite general viewpoint. If we look at a sys- 
tem of macroions and small ions (including counterions 
and salt ions) we can divide the contributions to the po- 
tential energy into three categories: Macroion-macroion 
(MM), macroion- ion (MI), and ion- ion (II). The rela- 
tive importance of these terms can be estimated based 
on simultaneous consideration of their magnitudes and 
distances, on which they are operating. Macroions re- 
pel each other, so that their interaction is of the or- 
der of (3umm ~ ^bZm exp (— nd)/d, the macroion-ion 
contribution (3umi ~ XbZmZi exp (—Ka)/a, and the 
ionic part f3uu rj XbZ] cxp (— K,di)/di. Here, Xb = 
e 2 / (AiresokBT) is the Bjerrum length, /3" 1 — ksT the in- 
verse of the thermal energy, ks the Boltzmann constant, 
T the temperature, Zm, Zi macroion and ion valence re- 
spectively, a the macroion radius, k = AnXsiZjCs+nZM) 
the screening parameter, d — n -1 / 3 and di = c s the 
mean macroion-macroion and ion-ion distance respec- 
tively, c s the salt content and n the macroion number 
density. Setting Z M = 1000, Zi = ±1, c sa it = 1 mM, 
a = 100 nm and the macroion volume fraction to 0.01, 
we get I3umm ~ 3 x 10~ 30 , (3u M i ~ 3 x 10 -4 , and 
/3uji « 3 x 10~ 2 . From this naive estimate, one can 
conclude that at na 3> 1 the two last contributions dom- 
inate the system thermodynamics. When the charge sign 
is taken into account, the negative MI contribution starts 
prevailing in the total energy, as the significant part of II 
contribution; consisting of nearly equal number of terms 
of opposite signs, cancels itself out. Due to strong screen- 
ing, the thermodynamic properties of such dispersion do 
not differ much from a simple electrolyte, except for small 
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layer of thickness k _1 around the macroion surface. One 
observes a much different picture in the regime of thick 
double layers ka C 1. A similar estimate for c sa i t = 1/iM 
gives (3umm ~ 80, 0umi ~ 5 and @un « 4 x I0~ 3 . 
If this energy per ion is weighted by the correspond- 
ing number of species, the MI contribution dominates 
so that the total Coulomb energy becomes negative [jj. 
We can also see that the II interactions are unlikely to 
influence the dispersion properties in case of monovalent 
ions. The first two contributions, however, have to be 
taken into account. The common method of dealing with 
this situation involves: (i) solving the Poisson-Boltzmann 
(PB) equation for a single macroion, (ii) renormalization 
of the MM interaction parameters based on the Debye- 
Hiickel-like approximation for the long distance part of 
the double layer, and (iii) solution of the one-component 
MM model with an effective interaction potential u e g 
[5 HT2I Il3| . This procedure respects the leading role of the 
MI interaction, while the effect of MM correlations in this 
approach enters only on the level of the one-component 
(OCM) description. The latter approach can be mod- 
ified in different ways to account for more pronounced 
role of macroion correlations and the resulting ch arge in- 
homogeneities using the Wigner-Seitz cell model [l3j. In 
this case, the structure of the double layer reflects the in- 
homogeneous macroion (and hence the counterion) dis- 
tribution via the cell construction. Numerical schemes 
based on the cell model and charge renormalization have 
been successful in describing properties of charged col- 
loidal dispersions [a, ll 2L Il3| . 

An alternative route to include macroionic contribu- 
tion into electrostatic screenin g is based on the jellium 
approximation for macroions lH El Although the 
range of validity of the jellium model might be limited 
to weakly correlated macroion systems, this model can 
be easily extended to different situations, i.e. rod-like 
colloids E3 or asymmetric electrolytes and, particularly, 
its equation of state takes a simple analytical form (see 
equation below). 

In this work, we follow the jellium description, which 
we improve using a simple construction that introduces 
spatial macroion correlations. The main goal of the 
present work is therefore to study the effect of macroion 
correlations on the parameters of the OCM and the 
corresponding dispersion thermodynamics. Another is- 
sue we would like to address is the calculation of the 
equation of state of the dispersion. It is well known 
that recovering the accurate equation of state basing on 
the OCM representation is problematic, as the ionic de- 
grees of freedom are omitted from the description and 
wall effects are not included on the macroion virial con- 
tribution |l8j . At the same time, the thermodynamic 
properties can be easily extracted from the same two- 
component (or multi-component) description that is used 
for calculation of the OCM effective charge and screening 
length, i.e. dire ctly from the WS cell or jellium model 
[ll El Q E El El H3. These models are usually 
solved using the non-linear PB equation including only 



one colloidal particle, or MC simulation 0, IT2I I2H |22| . 
In this work we will use non-linear PB equation with the 
jellium boundary conditions. 



II. MACROION CORRELATIONS ON THE 
MEAN-FIELD LEVEL 

Although the PB cell and the renormalized jellium 
models do not address the macroion degrees of freedom, 
they implicitly include a model of macroion distribution. 
The cell model supposes well separated particles, where 
the role of the neighboring macroions consists of limit- 
ing the volume available for small ions while jmmW 
is simply zero inside this cell. The double layer inside 
the cell is otherwise unperturbed by the rest of the sys- 
tem. In contrast, the jellium model assumes <7mm (?") = 1 
for r > 2a(diameter of the sphere), i.e. an ideal gas of 
macroions. As we already noted in the introduction, this 
approximation might be unsatisfactory for low-salt col- 
loidal dispersions, where macroions strongly repel each 
other even at the mean interparticle separation. We 
therefore will try to avoid considering uniform macroion 
distributions. On the simplest level, the uniform distri- 
bution can be replaced by a ^mm(^) taken in the form 
of a step function. This choice is motivated by the ob- 
servation that a charge-stabilized colloidal suspension at 
low salt shows a highly structured gMM(r) with a char- 
acteristic lengthscale described by the mean interparticle 
distance d = n~ 1 ^ 23]. In particular for r < d, gMAi(r) 
is almost zero, a feature that is known as the "correla- 
tion hole". The total charge density in the system at a 
distance r from the center of a macroion becomes 

p(r) = -Z BS eng M Ki{r) + ep + (r) - ep-(r), (1) 

where gMJu(r) = for r < d, and jmm^) — 1 for r > d; 
p±{r) are the concentrations of salt cations and anions 
and e is the elementary charge. We should stress that 
the "correlation hole" approximation for gMM(f) can be 
justified only for low-salt systems where the main peak 
position of gMM{r) scales with rt~ 1//3 . In a more gen- 
eral case, one should consider a " correlation hole" of size 
d* such that the main peak position scales with it and 
which should be valid for higher salt concentrations or 
weakly charged macroions. This point will be addressed 
elsewhere |24| . 

Our model includes now the major part of macroion 
correlations by placing a macroion at the center of its cor- 
relation hole with a size that depends explicitly on the 
concentration. The macroion distribution outside this 
hole is still approximated by an ideal gas of macroions, 
as in the jellium approach E3- We enforce the smeared- 
out background charge, Zb ac k, representing the other 
macroions around the tagged macroion to coincide with 
the effective charge, -Z e ff , as in the original work of Trizac 
and Levin The Zb ac k is determined by the elec- 

troneutrality condition for the total charge density in the 
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FIG. 1: Effect of the macroion correlation correction to the 
renormalized jellium model on the reduced pressure in a salt- 
free dispersion. The solid and dashed curves show the pres- 
sures in the moderate effective charge regime {ZmXb/cl = 10) 
predicted by the jellium and m-jellium, respectively and the 
solid and dashed curves with symbols the saturated effective 
charge regime {ZmXb/o, = 1000). The inset shows the rela- 
tive difference between the renormalized jellium and m-jellium 
pressures. This relative difference describes also the change 
in the effective charge in the salt-free system. 



bulk, 2c s sinh [e<p(oo) /ksT] = riZback, where </>(oo) rep- 
resents the electrostatic potential in the bulk. The clcc- 
troneutrality condition also allows us to determine the 
effective screening parameter, K 2 ff = Att Xb \J Z 2 s n 2 + 4c 2 
|l6j. We therefore call our model a modified jellium 
model (m-jellium). 

The effect of the introduced modification is illustrated 
in Fig. One can see that the pressure in a salt-free col- 
loidal dispersion becomes smaller in the m-jellium model 
as compared to the original jellium result at all macroion 
volume fractions. For both highly (Zj^Xb/o, — 1000) 
and moderately charged systems (ZmXb/ci = 10), the 
relative difference between the two models is maximal 
(about 0.3) at the lowest volume fraction and decreases 
as the volume fraction increases (it does not exceed 0.15 
at if = 0.1). The relatives differences between pressures 
from both models are illustrated in the inset of Fig. ^ 
The pressure in the salt-free case is proportional to the 
effective charge, equation (J2J, so the effective charge vari- 
ation in both models is also described by the same equa- 
tion. The difference between both models decreases with 
density because the size of the correlation hole also scales 



-1/3 



Therefore, it is expected that results from both 



as n 

models will coincide at higher volume fractions. 

For comparing different systems it is convenient to 
present the effective charge in the form Z^Xb/o,. In 
Fig. we compare the effective macroion charge as 

a function of its bare charge for the PB-cell, jellium and 
m-jellium models in the salt-free case. Each model pre- 
dicts different values of the effective charges at saturation 
(i?bare — > oo). However, for small bare charges all of them 
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FIG. 2: a) Effective charge for a salt-free suspension at a 
volume fraction of tp = 0.001 as a function of the bare charge 
and b) effective charges at saturation for charge-stabilized col- 
loidal suspensions, salt-free case, as a function of volume frac- 
tion. In both cases, dashed lines are for the PB-cell model, 
dotted lines for the jellium model and solid lines for the mod- 
ified jellium approach. 



recover the correct limiting behavior Z c s — Zbare- It is 
interesting to note that for the salt-free case at satura- 
tion the system properties are determined by only one 
parameter: the macroion volume fraction, ip = 47ra 3 n/3. 
In Fig. |2t>) we compare the saturated effective charges 
obtained from each model. The behavior of Z e ff in the 
range f < 10~ 2 can be understood by a compression of 
the ionic double layers. As if increases, counterions are 
pushed towards the macroion surfaces thus reflecting the 
win of the entropy over the energy and therefore leading 
to a gradual decrease of Z c ff until a well-defined minimum 
appears around (p « 10 -2 and we observe the growth in 
the effective charge at <p > 10~ 2 . This reduction in the 
effective charge in the m-jellium as compared to the orig- 
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After exclusion of the condensed counterions by the 
charge renormalization procedure only the free ions con- 
tribute to the pressure. In the salt-free case, c s = 0, 
the equation of state given by Eq. J2J simply reduces to 
(3P = n(l+Z e g), In the low electrostatic coupling regime 
(where Z e ff coincides with Zb are ), equation (J2J) recovers 
(see also Fig. EJ the ideal gas pressure f3P ~ p c (l + Zt, a re) 
[l3| . In Fig. we compare predictions from the renor- 
malized and the m-jellium models for the osmotic pres- 
sure data from Reus etal. [^J, obtained for deionized 
suspensions of bromopolystyrene particles. Although the 
overall agreement between experiments and both models 
is good, for ip < 0.07 the m-jellium approach performs 
visibly better while for higher densities the renormal- 
ized jellium model seems to be a better approximation. 
Further on, in Fig. ^p) we compare the results from 
both jellium-like models with data from primitive model 
simulations for salt-free asymmetric electrolyte with an 
asymmetry in charge 60 : 1 (we use the same cluster MC 
simulation protocol and settings as in Ref. Q with 80 
macroions) . Here we observe that both models reproduce 
the pressure behaviour qualitatively, while the m-jellium 
gives a better quantitative agreement in the wide range of 
macroion concentrations. An alternative representation 
of the same data is shown in the inset of Fig. JSb) , where 
both models are compared to the osmotic coefficients for 
the same simulated system. From this representation, it 
is clear that the m-jellium describes the results of nu- 
merical simulations better. Thus, our model improves 
the osmotic pressure for suspensions in the countcrion- 
dominated screening regime. 



FIG. 3: a) Pressure as a function of the volume fraction 
obtained from the m-jellium approach (solid line) and the 
renormalized jellium approximation (dashed line). Symbols 
correspond to the Reus' experiments |25| . No adjustable pa- 
rameters have been used and effective charges at saturation 
are considered, b) Pressure as a function of the volume frac- 
tion obtained from the m-jellium approach (solid line) and 
from the renormalized jellium approximation (dashed line) 
with Zbarc As/a = 19.47 and Xb/ci — 0.3245. Symbols corre- 
spond to primitive model simulations. 



inal jellium model follows from the weaker screening at 
r < rt~ 1//3 (now the macroions are excluded from the in- 
ternal parts of the double layer) , which leads to a higher 
free energy cost of charging the macroion. In other words, 
the weaker screening result to stronger attraction of the 
counterions to the macroion surface. 



IV. STRUCTURE 

The effective charge and the screening parameter com- 
puted from our m-jellium can be used to calculate the 
effective pair potential and the structure of the suspen- 
sion. We assume the effective pair interaction between 
macroions to have the Yukawa form [l3j 



(3u ef f{r) = Z 2 S X B 



exp(K cff a) 
1 + k cS a 



exp(-/t cff r) 



(3) 



The pair distribution of the macroions interacting 
through the effective pair potential @ can be computed 
using the Ornstein-Zernike (OZ) equation |26j. 



h MM (r) = c cff (r) + n / d\' c oS {r')h MM {\v - r'|), (4) 



III. PRESSURE 



Within the jellium model, the equation of state reads 



(]P = n+JZ 2 s n 2 +4c 2 s . 



(2) 



where hMM( r ) — 9mm(t) — 1 and c (r) is the so-called 
effective direct correlation function ^| [27]]. An addi- 
tional closure relation is needed to solve the OZ equa- 
tion. I n p articular, the Rogers- Young (RY) closure re- 
lation can be used to solve the OZ equation self- 
consistently. The RY closure enforces both thermody- 
namic and density fluctuations to be the same in order 
to calculate both the structure and thermodynamics of 
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a colloidal suspension and it is known to accurately de- 
scribe Yukawa systems [HI . Technically, this is done by 
computing the isothermal compressibility using the virial 
route for the OCM, Xy 1 = {df3PMM/dn) T , where Fmm 
is the macroion-macroion virial contribution and can be 
written as 

PPmm = Pc~ -r / g(r) rd r. (5) 

D Jr=2a ar 

The isothermal compressibility can be also computed 
through the relation Xc 1 = 1 — nc cS (q = 0) 0, where 
c oS (q) is the Fourier transform of the effective direct cor- 
relation function. Then, the RY closure relation enforces 
both routes to give the same isothermal compressibility 
(Xv = Xc) hi order to guarantee, at least partially, the 
thermodynamic consistency |28| . 

It is important to note that the OCM pressure is usu- 
ally very different from the pressure measured in the 
full multicomponent electrolyte. This discrepancy fol- 
lows from the dominance of microion contribution to the 
pressure. In fact, the total pressure can be well approx- 
imated by the partial pressure of the small ions mea- 
sured using the contact value theorem at the WS cell 
boundary or at infinity 0. Moreover, as we have seen in 
the last section, equation @ provides an excellent equa- 
tion of state for highly charged colloids. It is now clear 
that equation JSJ is not able to reproduce the isother- 
mal compressibility of the suspension and thus it cannot 
be used to reproduce the structure of the suspension. 
Nevertheless, if the total pressure (or compressibility) of 
the system is known from the original system (or a so- 
lution of the multicomponent m-jcllium model), then it 
can be combined with the OCM to get the thermody- 
namically consistent description. Thus, the OCM pres- 
sure should be replaced with the pressure in the full sys- 
tem, which is given by Eq. Q within our mean-field 
approach. We stress that the pressure in our procedure 
is evaluated from the same mean-field parameters that 
give the correct macroion structure. This remark con- 
cerns also the osmotic compressibility xt, which is re- 
lated to the infinite wavelength limit of the macroion 
structure factor |13| . Knowing the macroion structure 
is therefore enough to compute the correct osmotic com- 
pressibility of the multicomponent system. On the other 
hand, the isothermal osmotic compressibility can be cal- 
culated with Xv 1 = (9/3P/dn)T, where P is given by Eq. 
(0), with the same result. 

We would like to note also that these simple ideas 
can be applied for interpretation of experimental data. 
The effective macroion charge is usually extracted from 
the structure factors measured in scattering experiments 
through fitting the curve with a OZ-RY scheme with an 
effective Yukawa potential. A reverse version of our ap- 
proach can be used to predict accurately the equation of 
state of colloidal suspensions at low salt via the relation 
between the effective charge and the osmotic pressure. 

In figure 10^,) we compare the radial distribution func- 
tions (RDF) from our primitive model simulations (sym- 




FIG. 4: Macroion-macroion pair-correlation functions from 
primitive model simulation (symbols) and from our OZ-RY 
scheme (see text) with effective m-jellium parameters (solid 
lines), a) From the standard OZ-RY scheme (see text) with 
the OCM equation of state and PB-cell parameters (dashed 
lines), and b) from our OZ-RY scheme with effective jellium 
parameters (dotted lines). From left to right, the packing 
fractions are a) tp = 0.08, 0.01, 0.0025, 0.00125, and b) ip = 
0.089, 0.022, 0.0055, 0.0027, 0.00069. 



bols) with numerical results obtained using our OZ-RY 
scheme with the m-jellium parameters (solid lines). Also, 
results from the standard RY route |23| and screening 
parameters from the PB-cell model are provided (dashed 
lines). For the sake of clarity, the results with parame- 
ters from the renormalized jellium through the standard 
route are not shown, however it has found that they un- 
derestimate the structure (for a more detailed analysis 
see Ref. |29j). The highest macroion charge taken in 
simulations (Zbarc = 60) does not bring the system into 
fully saturated effective charge regime. The value for 
Xs/a, considered (among others) was 0.324. In figure 
Pfe,) we clearly observe a good agreement between sim- 
ulations and m-jellium results although still small differ- 
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ences around the main peak of the RDF for the higher 
volume fractions (ip > 0.04) can be observed. These small 
differences might result from the macroion shielding ef- 
fect, which usually rises the macroion distribution peak 
in dense suspensions |30j . However, our results show a 
visibly better agreement with simulations than those ob- 
tained with other mean-field approaches (dashed lines). 
In figure (0Jd) we compare the RDF from the primitive 
model simulations described in section ITTT1 with numeri- 
cal results obtained using our OZ-RY scheme with both 
m-jellium and jellium parameters. We observe that the 
jellium model always overestimates the structure of the 
suspension while m-jellium shows a better agreement. It 
is remarkable that the behavior of gum (r) in figure (jJJ is 
accurately reproduced by our scheme in the whole range 
of distances, whereas other mean-field schemes overes- 
timate the short-range behaviour of <7mm(^)- Also, we 
note that at high densities the mean peak of the pair 
correlation is predicted less accurately. 



The evaluation of colloidal effective charges in suspen- 
sions with weak screening is based on the solution to the 
nonlinear PB equation for the electric double layer in 
the presence of other macroions. Our model represents 
a modification of the renormalizcd jellium approach by 
Trizac and Levin. The obtained results suggest that these 
correlations are important in systems with thick double 
layers such as deionized colloidal suspensions, i.e. in the 
counterion-dominated screening regime. 

Our model describes well experimental results of ther- 
modynamic quantities in colloidal dispersions, such as 
the osmotic pressure. The static structure of colloidal 
suspensions is also accurately reproduced by an OCM 
scheme that uses the screening parameters of the m- 
jellium model. Moreover, from the experimental point of 
view, the effective parameters of the RY integral equation 
scheme can be used to predict accurately the equation of 
state via the relation between the effective charge and 
the osmotic pressure in our model. 



V. CONCLUSIONS 

We have introduced a procedure of including macroion 
correlations into a mean-field theory of screening in 
charged colloidal dispersions, which leads to modification 
of the effective parameters of the OCM: the macroion ef- 
fective charge and the screening length. 
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